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Abstract. Geometric phase phenomena in single neutrons have been observed 
in polarimeter and interferometer experiments. Interacting with static and time 
dependent magnetic fields, the state vectors acquire a geometric phase tied to the 
evolution within spin subspace. In a polarimeter experiment the non-additivity of 
quantum phases for mixed spin input states is observed. In a Si perfect-crystal 
interferometer experiment appearance of geometric phases, induced by interaction 
with an oscillating magnetic field, is verified. The total system is characterized by 
an entangled state, consisting of neutron and radiation fields, governed by a Jaynes- 
Cummings Hamiltonian. In addition, the influence of the geometric phase on a Bell 
measurement, expressed by the Clauser-Horne-Shimony-Holt (CHSH) inequality, is 
studied. It is demonstrated that the effect of geometric phase can be balanced by an 
appropriate change of Bell angles. 
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1. Introduction 

The total phase acquired during an evolution of a quantum system generally consists 
of two components: the usual dynamical phase (pd and the geometric phase <fi g . The 
dynamical phase, which depends on the dynamical properties, like energy or time is 
given by <pd = — V^/ H(t)dt. The peculiarity of the geometric phase lies in the fact 
that it does not depend on the dynamics of the system, but purely on the evolution path 
of the state. Considering a spin | system, the geometric phase is given by minus half 
the solid angle (<j) g = —Q/2) of the curve traced out. Since its discovery by M. V. Berry 
in 1984 pQ the topological concept has been widely expanded and has undergone several 
generalizations. 

The first experimental evidence of an adiabatic and cyclic geometric phase, 
commonly called Berry phase, was achieved with photons in 1986 [2] and later with 
neutrons [3]. Non-adiabatic [I] and non-cyclic [5] geometric phases as well as the off- 
diagonal case, where initial and final states are mutually orthogonal [6], have been 
considered. In addition to an early approach by Uhlmann [TJ, an alternative concept 
of geometric phase for mixed input states based on interferometry was developed by 
Sjoqvist et al. [8]. Here, each eigenvector of the initial density matrix independently 
acquires a geometric phase. The total mixed state phase is a weighted average of the 
individual phase factors. This concept is of great significance for all experimental 
situations or technical applications in which pure state theories oversimplify the 
description. Theoretical predictions have been tested using NMR and single-photon 
interferometry [5] [TO]. The idea has also been extended to the off-diagonal case [UJ [12]. 

Neutron interferometry [131 EH] provides a powerful tool for investigations of 
quantum phenomena. Particularly in the field of geometric phases, where the spatial 
as well as the spinor evolution leads to geometric phases: In the spatial case the two- 
dimensional Hilbert space is spanned by the two possible paths in the interferometer. 
It has been experimentally verified that a geometric phases for cyclic [15J, as well as 
non-cyclic evolutions [IB] , can be induced. In the case of spinor evolution, where the 
geometric phase is generated in spin subspace, the spinor rotations are carried out 
independently in each sub-beam, due to the macroscopic separation of the partial beams 
in the interferometer [T7J. Geometric phase effects are observed when the two sub-beams 
are recombined at the third plate of the interferometer followed by a spin analysis. For 
instance in [18], spin flippers in both beams clearly demarcate the separate contributions 
of the dynamical and geometric phase acquired in the spin subspace. 

The geometric phase in a single-particle system has been studied widely over the 
past two and a half decades. Nevertheless its effect on entangled quantum systems is 
less investigated. The Berry phase is an excellent candidate to be utilized for logic 
gate operations in quantum communication [15] . due to its robustness against noise. 
This has been tested recently using superconducting qubits |20j, and trapped polarized 
ultracold neutrons [21 J . Entanglement is the basis for quantum communication and 
quantum information processing. Therefore studies on systems combing both quantum 
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phenomena, the geometric phase and quantum entanglement, are of great importance 
[22| |23| |2l] . In the case of neutrons, entanglement is achieved between different degrees 
of freedom and not between different particles. Using neutron interferometry, with spin 
polarized neutrons, single-particle entanglement between the spinor and the spatial part 
of the neutron wave function |25j, as well as full tomographic state analysis [26], have 
already been accomplished. 

In this paper we report on miscellaneous geometric phase phenomena in neutron 
polarimetry as well as interferometry. In Section [2] polarimetric measurements of 
noncyclic geometric, dynamical and general phases are presented. In particular, our 
experiment demonstrates that the geometric and dynamical mixed state phases $ 9 
and resulting from separate measurements, are not additive [28] in the sense that 
the total phase resulting from a single, cumulative, measurement differs from $ 9 + <&d 
[291 [30]. Furthermore, we report on observation and manipulation of the geometric 
phase generated in one of the Hilbert spaces in a spin-path entangled single neutron 
system, namely the spin subspace. Section [3] focuses on the geometric phase generation, 
due to time-dependent interaction with a radio-frequency (rf) field. Here the system is 
characterized by an entangled state, consisting of neutron and radiation field, governed 
by a Jaynes-Cummings Hamiltonian. In Section H] the influence of the geometric phase 
on a spin-path entangled single neutron system is described. We demonstrate in detail 
how the geometric phase affects the Bell angle settings, required for a violation of a 
Bell-like inequality in the CHSH formalism. 

2. Experimental observation of non-additivity of mixed-state phases 

2.1. Neutron polarimeter scheme for phase measurement 

Consider the experimental setup shown in Fig. [TJ The polarizer P prepares the beam 
in the | ft) spin state. Subsequently, a coil carries out a 7r/2-rotation (Ui), creating 
a coherent superposition l/\/2(| ft) — i\ Jj.)) of spin eigenstates that acquire opposite 
dynamical phase due to Zeeman splitting within the field B z . Alternatively, one could 
say that the polarization vector r' rotates in the x, y plane after Ui. The second coil 
and some arbitrarily chosen propagation distance within B 2 implement a spin evolution 
for both eigenstates and thereby induces a pure state Pancharatnam (total) phase 4> 
[31]. The third coil (U{) carries out a — 7r/2-rotation in order to observe spin interference 
in the detector D after the analyzer A (both, P and A project the spin towards the 
+z direction). To obtain these interferences a phase shift rj is implemented by linear 
translation of the second coil. It was first stated in [32] that with such an apparatus 
one can obtain phases between spin eigenstates of neutrons, induced by a SU(2) 
transformation 

6, C) = e lS cos ^1 ft) <1r I - e~* C sin^l ft) (^ | + ^ sin f | (ft | + e~ iS cos £| ^ (^ | . (1) 

Equation (CQ) describes a general evolution of the system within static magnetic fields. 
The resulting total phase (p = arg(ft \U$\ ft) = 5 can be written as a function of the 
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Figure 1. a) Sketch of neutron polarimetry setup for phase measurement with overall 
guide field B z , polarizer P, three DC-coils to implement unitary operations U\, Tj\, 
U,!,, analyzer A and detector D. Greek letters denote polarization rotation angles. 
Shifting the second coil induces an additional dynamical phase ry/2 resulting in spin 
interference. Evolution of the state on the Bloch-sphere induced by 17$, associated 
to: b) Purely (noncyclic) geometric phase (2£ = ir/2). c) Combinations of dynamical 
and geometric phase on the Bloch sphere (0 < 2£ < 7r/2). 



maximum I max and minimum l m in intensity of the oscillations, exhibited by applying 
the phase shift rj. The intensity is proportional to cos 2 £ cos 2 8 + sin 2 £ cos 2 (£ — n). This 
only depends on the SU(2) parameters £, 8 and ( - set by choosing the spin rotation 
angles in the second coil and the additional propagation distance within the guide field 
B 2 , respectively. 

A neutron beam with incident purity r' = \ f'\ along the +z-axis (r' = (0,0, r')) 
is described by the density operator p» n (r) = 1/2(1 + r'a z ). For mixed input states, 
< r' < 1. In this case [33] we find the intensity oscillations to be proportional to 
1 — r' 

P = — — + r' (cos 2 £ cos 2 8 + sin 2 £ cos 2 (C -rj))- ( 2 ) 

Considering again the maxima and minima of the intensity oscillations, one obtains the 
mixed state phase 

*/ i\ _ I iCiJln - 1/2(1 - r')]/r' 

nr ) -ccos y + ^ _ + _ ^ _ ^ { 6) 

with a normalization factor I? = 21^/(1 +r'). Iq is the intensity measured with = 1. 

Generally, the noncyclic geometric phase is given by <p g = —Q/2, where Q is the 
solid angle enclosed by an evolution path and its shortest geodesic closure on the Bloch 
sphere j5]: <f> g and the total phase are related to the path by the polar and azimuthal 
angles 2£ and 28 respectively, so that the pure state geometric phase can be written as 

s = 0-0 d = 5[l-cos(2O]. (4) 
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Figure 2. Mixed state phases determined from measured intensity oscillations using 
Eq. ([5]). Dynamical phase (a), geometric phase (b) versus input purity r for three 
evolution paths, i. e. three settings of the second coil (angle 2£) and flight distance 
after it (angle 25). The legends indicate evolutions. Solid lines are theory curves using 
the rightmost data points as reference. 



4>d is the dynamical phase. By proper choice of 2£ and 25, can be set to generate 
purely geometric, purely dynamical, or arbitrary combinations of both phases, e.g. in 
Figs. [lb and c: For instance, we can choose to induce a purely geometric phase by 
selecting 2£ to be equal to tt/2. 

The theoretical prediction for the mixed state phase is [8j [33] 

$ = arctan (r' tan 5) (5) 

Note that Eq. (jSJ only depends on the parameter 5 and the purity r' . Again, as can be 
seen also from Eq. @, the parameter £ only determines the portion of dynamical phase 
4>d contained in the total phase 0. 



2.2. Experiment 

To access Eq. flH]) experimentally r' has to be varied. In addition to the DC current, 
which effects the transformation Ui, random noise is applied to the first coil, thereby 
changing B x in time. Neutrons, which are part of the ensemble Pi n (r'), arrive at different 
times at the coil and experience different magnetic field strengths. We are left with the 
system in a mixed state r = (0, — r, 0) where r < 1 |34j. 

A neutron beam - incident from a pyrolytic graphite crystal - with a mean 
wavelength A ~ 1.98 A and spectral width AA/A ~ 0.015, was polarized up to 99% 
by reflection from a bent Co-Ti supermirror array. The final maximum intensity was 
about 150counts/s at a beam cross-section of roughly 1cm 2 . A 3 He gas detector was 
used. Spin rotations around the +x-axis were implemented by magnetic fields B x from 
DC coils on frames with rectangular profile (7x7x2 cm 3 ). B z was realized by two 
rectangular coils of 150 cm length in Helmholtz geometry. Low coil currents of about 
2 A corresponding to field strengths of up to 1 mT were required for the spin rotations 
and to prevent unwanted depolarization. The noise from a standard signal generator 
consisted of random DC offsets varying at a rate of 20 kHz. The experimental data 
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Figure 3. Filled markers: Measured total mixed state phase <£>t f versus purity r 
for two examples of XJtot associated to the total pure state phases <f> g + cf>d (see text). 
Open markers: <£> g + &d as calculated from measured data. Filled (empty) bars show 
measured mixed-state geometric (dynamical) phases. The solid and dotted theory 
curves assume cither non-additivity or additivity, respectively. 



shown in Fig. |2] reproduce well the r'-dependence predicted by Eq. flH]). 
2.3. Non-additivity 

Our experiment focuses on a special property of the mixed state phase: its non- 
additivity. The Sjoqvist mixed state phase [8] is defined as a weighted average of 
phase factors rather than one of phases. So it is true only for pure states that phases 
accumulated in separate experiments can be added up to the usual total phase in the 
following sense: Suppose a geometric pure state phase 4> g is induced in a first, and a 
dynamic pure state phase 4>d in a second measurement. Applying (pEJ) we can also choose 
a combination of angles 2£ and 26 leading to a transformation U iot , so that we measure 
the total pure state phase 4> g + 4>d in a third experiment. However, the result of this latter 
experiment for the system initially in a mixed state is &tot{r) = arctan [rtan(</> 5 + <pd)}- 
The total phase is not equal to + $d(r), with $ 9 (r) = arctan (r tan <p g ) and 

< ^ ) d( r ) — arctan (r tan 0^). Two examples of data confirming this prediction are shown 
in Fig. |3j 

3. Geometric phase generation in an oscillating magnetic field 

The evolution of a system consisting of neutron, static magnetic field and quantized 
rf-field can be described by a photon-neutron state vector, which is an eigenvector of 
a Jaynes-Cummings (J-C) Hamiltonian [35l [36], adopted for this particular physical 
configuration [37]. Since two rf- fields (the reason for the second rf-field is explained in 
Section I3TT]) . operating at frequencies u and oj/2, are involved in the actual experiment, 
the modified corresponding J-C Hamiltonian is denoted as 

%j-c = -p— V 2 - iiB (r)a z + h^a^a^ + — a^a^) (6) 



New Aspects of Geometric Phases in Experiments with polarized Neutrons 7 

(R^) ( \ Ej( a, /2)/ \ \ 

with a = \(o~ x + io~ y ). The first term accounts for the kinetic energy of the neutron. 
The second term leads to the usual Zeeman splitting of 2\h\Bq. The third term adds the 
photon energy of the oscillating fields of frequencies uj and cu/2, by use of the creation 
and annihilation operators a) and a. Finally, the last term represents the coupling 
between photons and the neutron, where N Uj = (o^ a Uj ) represents the mean number 
of photons with frequencies ojj in the rf-field. 

The state vectors of the oscillating fields are represented by coherent states \a), 
which are eigenstates of and a. The eigenvalues of coherent states are complex 
numbers, so one can write a\a) = a\a) = |a:|e 49i |a)with|a:| = yN. Neutrons interacting 
with electromagnetic quanta are usually described by the 'dressed-particle' formalism 
[37] . in analogy to the dressed-atom concept [38J developed nearly two decades before. 
Using Eq. ((7j) one can define a total state vector including not only the neutron system 
I^n)) but also the two quantized oscillating magnetic fields: 

l^tot) = K) <g> \a u/2 ) ® |^n)- (8) 

In a perfect Si-crystal neutron interferometer the wavefunction behind the first plate, 
acting as a beam splitter, is a linear superposition of the sub-beams belonging to the 
right (| I)) and the left path (|II)), which are laterally separated by several centimeters. 
The sub-beams are recombined at the third crystal plate and the wave function in 
forward direction then reads as |^n) ^ |I) + |H), where |I) and |II) only differ by an 
adjustable phase factor e %x (x = —N ps b c XD, with the atom number density iVp S in the 
phase shifter plate, the coherent scattering length b c , the neutron wavelength A and the 
thickness of the phase shifter plate D). By rotating the plate, \ can be varied. This 
yields the well known sinusoidal intensity oscillations of the two beams emerging behind 
the interferometer, usually denoted as O- and H-beam [T3] . 

In our experiment, only the beam in path II is exposed to the rf-field of frequency 
oj, resulting in a spin flip (see Fig.|U(a)). Interacting with a time-dependent magnetic 
field, the total energy of the neutron is no longer conserved after the spin-flip 
[39| HOl SH H21 H31 Hi]. Photons of energy fvjj are exchanged with the rf-field. 

The time-dependent entangled state, which emerges from a coherent superposition 
of 1 1) and | II), is expressed as 

|*tot) = K) <B> K /2 ) ® ^=(|I) + e^e^|II) ® e<*\ ^)), (9) 

for a more detailed description of the generation of |^ t ot) see [^5] - 

The effect of the first rf-flipper, placed inside the interferometer (path II), is 
described by the unitary operator U(<pi), which induces a spinor rotation from | ff) 
to | Jj), we denoted U((f>i)\ ft) = e^ : | JJ-). The rotation axis encloses an angle <pi with 
the x-direction, in the rotating frame, and is determined by the oscillating magnetic 
field flW = cos(ojt + 0i) • y. Formally one can insert a unity operator, given by 
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Figure 4. (a) The experimental apparatus for observation of geometric phase. The 
spin state acquires a geometric phase 7 during the interaction with the two rf-ficlds 
and is flipped twice. Finally, the spin is rotated by an angle 8 = n/2 (in the x,z 
plane), by a dc-spin turner, for a polarization analysis and count rate detection, (b) 
The Bloch-sphere description depicts the acquired geometric phase given by minus half 
the solid angle depending on the phase <j)\ of the rf-field. The effect of the beam block 
is explained in Section 2) 



]l = m(0o)Wo), yielding 

u(<k)\ fr) = u(<k) ^MUM \ ft) = e n l (10) 

where U{<fto) can be interpreted as a rotation from | "f|~) to | JJ-), with the x-direction 
being the rotation axis (0o = 0), and t/f(0 o ) describes a rotation about the same axis 
back to the initial state | ff). Consequently, U{4>i) 4>o) can be identified to induce the 
geometric phase 7, along the reversed evolution path characterized by 0o (| J|) to | 
followed by another path determined by (pi (| ft) to | JJ.)), see Fig.|H(b). In the rotating 
frame of reference [16] the two semi-great circles enclose an angle <pi and the solid angle 
VL = —20:, yielding a pure geometric phase 

7 = "^/2 = 0i, (11) 

which is depicted in Fig.H](b). The entangled state, as described in [22J, is represented 
by 

IW7)> = 7f (l 1 ) ® 1 + 1 11 ) ® en l ^>). ( 12 ) 

including the geometric phase 7 = <p\. 



3.1. Experimental Setup 

As in a previous experiment [45] . the spin in path |II) is flipped by a rf- flipper, which 
requires two magnetic fields: A static field B ■ z and a perpendicular oscillating field 
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Figure 5. (a) Typical interference patterns of the O-beam, when rotating the phase 
shifter plate (x)-(b) A phase shift occurs by varying <j>i determining the geometric phase 
7 . The sign of the geometric phase 7 ± depends on the chosen initial polarization. 



B^ = Bj[' cos(ut + 0i) • y satisfying the amplitude and frequency resonance condition 

n(«) 2\fJ,\B B\ 

B - =7H andW = ^ (1+ 165f ) ' (13) 

where /i is the magnetic moment of the neutron and r denotes the time the neutron is 
exposed to the rf-field. The second term in cu is due to the Bloch-Siegert shift [17]. The 
oscillating field is produced by a water-cooled rf-coil with a length of 2 cm, operating at 
a frequency of u/2n = 58 kHz. The static field is provided by a uniform magnetic guide 
field B^ ~ 2 mT, produced by a pair of water-cooled Helmholtz coils. 

The O-beam passes the second rf-flipper, operating at oj/2tt = 29 kHz, which 
is half the frequency of the first rf-flipper. The oscillating field is denoted as 
B^ 2 ' ) cos ((oj/2)t+(f)ii) -y, and the strength of the guide field was tuned to B^^ ~ 1 mT 
in order to satisfy the frequency resonance condition. This flipper compensates the 
energy difference between the components from the two interfering paths, by absorption 
and emission of photons of energy E = hu/2. By choosing a frequency of oj/2 for 
the second rf-flipper, the time-dependence of the state vector is eliminated since both 
components acquire a phase e ±Ja; / 2 (*+ T ) ( depending on the spin orientation. Only a 
constant phase offset of e^ 7 , where T is the propagation time between the centre of 
the first and second flipper coil, remains in the stationary state vector. This phase 
contribution, together with a dynamical phase contribution, resulting from Larmor 
precession within the guide field regions B^ and B^^ 2 " 1 (pointing in + z-direction), 
are omitted here because they remain constant during the entire experiment. Finally, 
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the spin is rotated by an angle S = it/ 2 (in the x, z plane) with a static field spin- 
turner, and analyzed due to the spin dependent reflection within a Co-Ti multi-layer 
supermirror along the z-direction. Intensity oscillations in forward direction (O-beam) 
are plotted in Fig. 0(a). 

In a non- dispersive arrangement of the monochromator and the interferometer 
crystal the angular separation can be used such that only the spin-up (or spin- 
down) component fulfils the Bragg-condition at the first interferometer plate (beam 
splitter). Therefore it is possible to invert the initial polarization simply by rotating the 
interferometer by a few seconds of arc, which is expected to lead to an inversion of the 
geometric phase. Figure [5(b) shows a plot of the geometric phase A7 =b versus <fii, with 
0n = 0. As expected, the slope s is positive for initial spin up orientation (s = 1.007(8)), 
and negative for the spin down case (s = —0.997(5)), as predicted in Eq.tfTTT) 

4. Geometric phase effects on a spin-path entangled System 

In this Section the influence of the geometric phase on a Bell measurement [IB] , expressed 
by the Clauser-Horne-Shimony-Holt (CHSH) [19] inequality, as proposed in [22] is 
discussed . Following the notation given in [22], the neutron's wavefunction is defined 
via tensor product of two Hilbert spaces: One Hilbert space is spanned by two possible 
paths in the interferometer given by |I) and |II); the other one by spin-up and spin-down 
eigenstates, denoted as | ff) and | Jj), with respect to a quantization axis along a static 
magnetic field. For this experiment we focus on the neutron part of Eq.([9]) and omit all 
phases but the geometric phase 7: 

|* N ( 7 )> = -J=(|I) <g> I fr) + |II) ® e*| ^>). (14) 

As in common Bell experiments a joint measurement for spin and path is performed, 
thereby applying the projection operators for the path 

Pl(a) = \±a)(±a\, (15) 

with 

I +a) = cosy|I) +e ia2 siny|II) and | - a) = - sin y |I) + e iQ2 cos y |II), (16) 

where a± denotes the polar angle and «2 the azimuthal angle for the path. This is done 
in analogous manner for the spin subspace with ft as the polar angle and ft as the 
azimuthal angle for the spin. Introducing the observables 

A p (a) = P$(a) - P p (a) and B s (/3) = P s + (/3) - P b l(/3) (17) 

one can define an expectation value for a joint measurement of spin and path along the 
directions a and (3 

E(a,(3) = (ty\A p (a) g> B s (f3)\^} = -cos«icos/3i - cos(a 2 - fh + 7) sin a x sin ft (18) 
= - cos(cti - ft) for (ot 2 - fa) = —7- 
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Next, a Bell-like inequality in CHSH- formalism [49J is introduced, consisting of four 
expectation values with the associated directions ex, ct' and (3, for joint measurements 
of spin and path, respectively 

S(a, cx', (3, (3', 7 ) = \E(cx, (3) - E(a, (3') + E{d , (3) + E(cx', (3') | (19) 

Without loss of generality one angle can be eliminated by setting, e.g., a = 
(«i = a 2 = 0), which gives 

S(a', (3, /3' , 7) = — sin o>[ ^cos(o4 — Pi — 7) sin [3\ + cos(a 2 — P' 2 — 7) sin f3[) (20) 
— cos a^(cos (3i + cos f3[) — cos j3\ + cos (3[ . 



The boundary of Eq. (1191) is given by the value 2 for any noncontextual hidden- variable 
theories [50]. Keeping the polar angles a[, (3i and (3[ constant at the usual Bell angles 
a[ = |, Pi — |, /3( = ^ (and azimuthal parts fixed at a' 2 = (3 2 = (3' 2 = 0) reduces S to 

5(7) = I - V2 - V2 cos 7 |, (21) 

where the familiar maximum value of 2 a/2 is reached for 7 = 0. For 7 = n the value of 
S approaches zero. 

4-1. Polar Angle Adjustment 

Here we consider the case when the azimuthal angles are kept constant, e.g., a' 2 = f3 2 = 
(3' 2 = (a 2 = 0), which is depicted in Fig. [6] The S'-function reads as 

S(a[, Pi, f3[, 7) = — sin a'i (^ cos 7 srn Pi + cos 7 srn P'xj ~ cos (cos /?i + cos P[) (22) 

— COS /?i + cos P[ . 

The polar Bell angles Pi, P[ and a[ (ai = 0), yielding a maximum S- value, can be 
determined, with respect to the geometric phase 7, by calculating the partial derivatives 



PATHoc 2 =a 2 '=0: ^(a^O) SPIN ^=^2 = 0:. ITTXTTl 



x =0 (a^^/2) 




= jt (o / 1 - L =3n/2) 



IIIXIII (a^- = 7t) 




Pi" (y) 



Pf (y) 



Figure 6. Bloch-sphere description includes the measurement settings of a and /3(<5), 
determining the projection operators, used for joint measurement of spin and path. 
a is tuned by a combination of the phase shifter (x) an d the beam block, and (3 is 
adjusted by the angle 5. 
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Figure 7. (a) Polar adjusted ^-values versus geometric phase 7 with adapted Bell 
angles (0% and j3[) according to the geometric phase 7. (b) the corresponding modified 
Bell angles are plotted versus the geometric phase 7. 



(the extremum condition) of S in Eq.( l22p (see [22] for more elaborated deduction): The 
solutions are given by 

7T 

/3\ = arctan (cos 7), f3[ = it — (5\ and ot x = — , (23) 

which are plotted in Fig.[7](b) (denoted as theoretical predictions). With these angles 
the maximal S decreases from S = 2y/2 for 7 : — > | and touches at 7 = ~ even the 
limit of the CHSH inequality S — 2. Within the interval 7 : ~ — > ir the value of S 
increases again and returns to the familiar value S = 2\/2 at 7 = it. 

Experimentally, the probabilities of joint (projective) measurements are propor- 
tional to the following count rates iV^ with (i,j = +, — ), detected after path (a) and 
spin (j3) manipulation: 

F( m = iV ++ K g) - iV + _K (3) - N_ + ((x, (3) + M_(q, g) 
l "' Pj N ++ ( a ,f3) + N + 4cx,f3) + N^( a ,f3) + N^( a ,f3y { ] 
with for example 

iV ++ (a,/3) = N ++ (cx, (A,0)) oc (* N ( 7 )|pP(a) ® P^(/?i, 0)|* N ( T )) (25) 

In the case of iV + _(a, /3) the count rate is given by A^ ++ (a, (/Jj 1 , 0)), where (3± = (3i + n. 

The procedure is applied to the count rates N + -(a, (3) and iV (ct,/3). With these 

expectation values S can be calculated as defined in Eq. ([T9]) . 

Projective measurements are performed on parallel planes defined by 02 = a' 2 = 
P2 — P'2 = 0- For the path measurement the directions are given by a : 07 = 0, cti = , 
and a' : a[ = ir/2, a' 2 = 0. 

The angle a, which corresponds to + z (and — z for = a\ + n = n, a 2 = 0) is 
achieved by the use of a beam block which is inserted to stop beam II (I) in order to 
measure along +z (and — z). The corresponding operators are given by 

Pt(«i = 0,a 2 = 0) =|I)(I| 

P p z (ai = n,a 2 = 0) = |II)(II|, (26) 
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The angle a! is set by a superposition of equal portions of |I) and |II), represented 
on the equator of the Bloch sphere. The interferograms are achieved by a rotation of 
the phase shifter plate, associated with a variation of the path phase x- All path scans 
are repeated at different values of the spin analysis direction S in order to determine 
Pi and (3[ for a maximal violation of the Bell-like CHSH inequality. The projective 
measurement for a[ = n/2, a' 2 = corresponds to a phase shifter position of x = (and 
a[ ±= a[ + Ti = 37r/2, a 2 = to x — n )- Projection operators read as 

PIM = |>«2 = 0) = + |H»«I| + (H|)) (27) 

P^M 1 - = y >4 = °) = ~ ~ (28) 

Using the measurement curves from Eq. (125|) and Eq. (|2T|) . the S- value is calculated 
according to Eq. ffT9]) as a function of the parameters (3i and (3[, which are varied 
independently. The local maximum of the S(P[,Pi) is determined numerically and 
plotted in Fig.[7](a), with the corresponding values for f3\ and (3[ in Fig.d(b). 

4-2. Azimuthal Angle Adjustment 

Next we discuss the situation where the standard maximal value S = 2^/2 can be 
achieved by keeping the polar angles a[, fi\ and (3[ constant at the Bell angles a[ = f , 
f3\ = f , /3[ — ^jr, (a i = 0), while the azimuthal parts, a' 2 , P>2 and (3' 2 (a 2 = 0), are 
varied. A Bloch sphere description of this configuration can be seen in Fig. [HI The 
corresponding S function is denoted as 

S(a' 2 , (3 2 , (3> 2 , ,) = \-V2-f (cos(o4 - fa - 7 ) + cos(^ - (3' 2 - 7 )) | . (29) 

The maximum value 2 a/2 is reached for 

f3 2 = 02, and a 2 — (3' 2 = 7 (mod7r). (30) 

For convenience (3 2 = is chosen. 

Experimentally the angle between the measurement planes is adjusted by one 
azimuthal angle (a' 2 ), which is deduced by phase shifter (%) scans. 



PATHcc 2 = 0: f |I><l| (a, = 0) 




5PIN p 2 =p2 = 0: 



|n><n| (a^ = 7t) pi=3n/4 




[txtri 

p'^ = - Tt/4 



Figure 8. Bloch-sphere description includes the measurement settings of a and (3(6), 
determining the projection operators, used for joint measurement of spin and path. 
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Figure 9. (Color online) (a) Azimuthal adjusted S'-values versus geometric phase 
7 with balanced Bell angle (a' 2 ) according to the geometric phase 7, and without 
corrections, (b) the corresponding modified Bell angle is plotted versus the geometric 
phase 7. 



For the spin measurement the directions are fixed and given by (3: $\ = 7r/4, (3 2 = 
and (3': (3[ = 3vr/4, f3' 2 = (together with ^ = -3vr/4, f}[ x = -it /A. For the projective 
path measurement the fixed directions read as a% = (aj- = tt), for measurements with 
beam block, and a[ = tt/2 (a' L = 3tt/2). Phase shifter (x) scans are performed in 
order to determine a' 2 for a maximal violation of the Bell-like CHSH inequality yielding 

S = 2y/2. 

As predicted by Eq. (l3"0i) the constant maximal S value of 2\/2 (see Fig.|9](a)) is 
found for a' 2 = 7, which is plotted in Fig.|9](b). In Fig.[9](a), the case is included 
where no corrections are applied to the Bell angles. According to Eq.( |2TT) the familiar 
maximum value of 2\/2 is reached only for 7 = 0, and at 7 = it the value of S = is 
found. 

This experiment demonstrates in particular, that a geometric phase in one subspace 
does not lead to a loss of entanglement. Two schemes, polar and azimuthal adjustment 
of the Bell angles, are realized, balancing the influence of the geometric phase. The 
former scheme yields a sinusoidal oscillation of the correlation function S, such that 
it varies in the range between 2 and 2\/2 and, therefore, always exceeds the boundary 
value 2 between quantum mechanical and noncontextual hidden-variable theories. The 
latter scheme results in a constant, maximal violation of the Bell-like CHSH inequality, 
where S remains 2a/2 independent of the value of the geometric phase 7. 

5. Concluding remarks 

Neutron optical experiments are used for studying characteristics of phases of geometric 
origin. First, non-additivity of the mixed state phase has been observed in a polarimetric 
experiment. Since the purity of quantum states in real experiments is always smaller 
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than 1, non-additivity is of importance in all applications of quantum phases. Thinking 
about phase gates, it means that the purity of the utilized quantum system has to 
be considered when inducing phases for further processing. Second, a technique for 
geometric phase generation has been established by means of a precise spin manipulation 
due to interaction with rf-fields, in an interferometric setup. Applying the formalism of 
the Jaynes-Cumming Hamiltonian to the patterns in the observed outgoing beam of the 
interferometer, we find good agreement between experiment and theory. This technique 
is also applied to phase manipulations of the spin subspace in a triple-entanglement 
experiment with neutrons, which will be topic of a forthcoming publication. Finally, 
the effect of the geometric phase on the entanglement of the system, has been analyzed 
in detail, using a Bell-like CHSH inequality. It is demonstrated, how the effects of 
the geometric phase on the outcome of a Bell measurement can be balanced by an 
appropriate change of Bell angles. Neutrons have proved to be a suitable quantum 
system for studying topological effects. Interferometric as well as polarimetric techniques 
will lead to further investigations, relevant for possible applications of the geometric 
phase. For instance, we are planing a polarimetric experiment, in which the geometric 
phase for non-unitary evolutions is considered. 
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